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MODELS FOR THE UV-X EMISSION IN AGN: INFLUENCE OF 
DIFFERENT PARAMETERS 
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We present the results of computations performed with a new photoionisation-transfer code designed for hot 
Compton thick media. Presently the code solves the transfer of the continuum with the Accelerated Lambda 
Iteration method (ALI) and that of the lines in a two stream Eddington approximation, without using the local 
escape probability formalism to approximate the line transfer. We show the influence that the approximations 
made in the treatment of the transfer and of the atomic data can have on the broad band UV-X spectrum and 
on the detailed spectral features (X-ray lines and ionization edges) emitted/absorbed by an X-ray irradiated 
medium. This transfer code is coupled with a Monte Carlo code which allows to take into account direct and 
inverse Compton diffusions, and to compute the spectrum emitted up to hundreds of keV energies. The influence 
of a few physical parameters is shown, and the importance of the density and pressure distribution (constant 
density, pressure equilibrium) is discussed. 



1. The codes 

Using a photoionisation-transfer code specially 
designed for hot Compton thick media (T > 
W^K, Thomson depth Tes can be as large as a 
few hundreds), we are able to compute the ther- 
mal and ionisation structure and the spectrum 
emitted and reflected by a slab of gas illuminated 
on one side or on both sides by a given spectrum 
in a plane-parallel geometry (for a description of 
the code, see Dumont et al. [ Q and Dumont & 
Collin [ H). Contrary to the other photoion- 
isation codes which use the escape probability 
(EP) formalism, our code solves the transfer of 
the continuum with the Accelerated Lambda It- 
eration (ALI) procedure and that of the lines with 
a two stream approximation (ALI is being imple- 
mented presently for the lines). It is coupled with 
a Monte-Carlo code which allows to take into ac- 
count Compton and inverse Compton diffusions, 
and to compute the spectrum emitted up to hun- 
dreds of keV energies, in any geometry. 

The gas composition includes 10 elements and 
102 ions. H is treated as a 6-level atom, H-likc 
and Li-like as 5-level atoms. He-like ions as 8- 
levels atoms and matrix inversions are performed 
to compute the levels populations. The other 
ions are treated roughly as two level atoms plus 
a continuum. All coUisional and radiative pro- 



cesses (including ionizations from and recombina- 
tions onto excited levels) are taken into account. 
Recombinations onto upper levels (n > 6 for H, 
n > 5 for H- and Li-like ions, n > 8 for He-like 
ions) are added to our highest excited level. 

2. Results 

We consider here two cases : 

• the density is constant across the slab 

• the pressure is constant across the slab 

In the present examples the incident continuum 
is a power law with a slope of unity (in flux, F,^) 
from 0.1 eV to 100 keV. We take cosmic abun- 
dances. The density at the surface of the slab is 
lO^^cm^'^ and the incident flux is equal to 8 10^'^ 
erg cm""^ s^^ (corresponding to an ionisation pa- 
rameter ^ — L/{nHR^) — 1000 erg. cm. s^^) in 
constant pressure and constant density cases. 

2.1. Influence of the transfer treatment 

In Dumont et al. [ Q, it was shown that the 
neglect of the returning radiation often used in 
photoionisation codes, i.e. the "outward only" 
approximation, leads to strong temperature er- 
rors, in particular at the illuminated face of the 
slab when the slab is thick {t^s > 0.1). 
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Figure 1. Radiative rate cooling- heating curves 
(full lines) and net cooling curves due to lines 
(dashed lines) for different optical thicknesses in 
the constant pressure case (t^s — and equilib- 
rium temperature T^q — 2 lO^K, t^s — 0.64 and 
Teq = 3.5 10'* if). One can see that the line cool- 
ing which dominates the energy balance becomes 
negative close to the equilibrium temperature. 



One other difference between EP and transfer 
lies in the fact that the line emissivity can be 
negative whereas it is impossible in EP (see [ 
|l[). This translates into a difference in the energy 
balance, which can have important consequences 
on the equilibrium temperature (see Fig. |^ for an 
example in the constant pressure case). 

2.2. Influence of the atomic data 

The model atoms are important not only for 
the thermal balance, but also for the ionization 
equilibrium. Fig. ^ shows the temperature profile 
and Fig. ^ the fractional abundances of oxygen 
versus optical thickness, for two different treat- 
ments of the atomic data, in the constant den- 
sity case. When only H-like ions are treated with 
multi-level atoms, the temperature is higher than 
in the treatment with H-like, He-like and Li-like 
multi-level ions for an optical thickness smaller 
than about 2 and then is lower, for larger optical 
thicknesses. 

Treating H-, He-, Li- like ions as multi-level 
atoms produce more lines than with only H-like 
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Figure 2. Temperature profile versus optical 
thickness, for a treatment with only multi-level 
H-like ions (dashed lines), and with H-like, He- 
like and Li-likc multi-level ions (full lines) for a 
constant density model. 



multi-level atoms; all these added lines cool the 
medium which leads to a cooler temperature in 
the case of H-, He-, Li-like multi-level ions. The 
temperature profile declines faster in the H-like 
description than in the 'total' description between 
Tes = 2 and t^s = 3 because in that place, pho- 
toionization from the excited levels of Li- like ions 
are as important as those from the ground level in 
contradistinction to H-like and He-like ions where 
the photoionisation from excited levels is always 
negligible. When looking at the fractional abun- 
dances of oxygen, one sees that O VHI is present 
deeper in the cloud in the description with H-, 
He-, Li-like multi level atoms than in the treat- 
ment with only H-like multi level atoms. This is 
due to an approximate treatment of losses due to 
recombination radiation in the two-level atoms. 

2.3. Comparison of constant density and 
constant pressure models 
The slab being illuminated with the same irra- 
diating flux, we compare the results in the con- 
stant density and in the constant pressure cases. 
The temperature profiles and fractional abun- 
dances of oxygen are shown on Figs. |^, ^ and 
||, H respectively. At first sight, one sees that the 
thermal and ionization structure are completely 
different in the two different cases, owing to the 
phase change of the photoionized gas in the con- 
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Figure 3. Fractional abundances of oxygen versus 
optical thickness for the same two treatments of 
the atomic data. As in Fig. ||, the example shown 
here is for a constant density model. 



Figure 5. Fractional abundances of oxygen versus 
optical thickness for the constant pressure case. 
Note that, in the constant density case (fig 
only highly ionized species are present in the re- 
gions contributing to the emission, while in the 
constant pressure case, both high and low ioniza- 
tion species are present. 




Figure 4. Temperature profile versus optical 
thickness for the constant density and the con- 
stant pressure cases. The temperature decreases 
abruptly in the constant pressure case, when the 
soft X-ray photons are absorbed. 
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Figure 6. Comparison of the reflected spectrum 
between 1 eV and 25 keV in a constant density 
and a constant pressure cases. 
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We intend to treat all ions as multi-level atoms 




and to add levels to those already treated like 
this. We will take into account some forbidden 
lines which could be important for the cooling of 
the medium in some places. 

Though there are still many improvements to 
perform, we already use theses codes to describe 
AGN spectra in the UV and X-ray range, in their 
validity range. In particular, the code is used to 
compute the spectrum emitted by an irradiated 
slab in hydrostatic equilibrium (see Rozahska et 
al. in this workshop). 
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Figure 7. Zoom of Fig. ^ between 100 eV and 
20 keV. A few lines and edges are labelled in the 
figure, but others are identified in Table 111 



stant pressure case (see Krolik et al. [ ^) when 
soft X-ray photons have been absorbed. The tem- 
perature decreases abruptly and consequently the 
density increases in the constant presssure case. 
Accordingly, the emission-reflection spectra are 
different. Fig. ^ shows the reflected spectra in the 
constant density and the constant pressure cases, 
and Fig. |^ is a zoom of Fig. ^ between 0.1 and 20 
keV on which a few lines and edges are identified. 
Some line equivalent widths (EWs) of interest, in 
the constant density and constant pressure cases, 
are given m Table |l|. We recall that the EWs 
are those measured only on the reflected spectra 
with no dilution due to the incident spectrum. 

Due to the absence of layers with an interme- 
diate temperature in the constant pressure case, 
the "intermediate" ions are almost absent. As 
a consequence, the iron lines correspond to both 
high and low ionization species, in contradistinc- 
tion to the constant density case. Finally, due to 
the absence of intermediate temperature layers 
giving rise to soft X-ray reflection, the reflected 
spectrum in the constant pressure case displays a 
narrower "Blue Bump" than in the constant den- 
sity case. 

3. Foreseen improvements 

Presently we are implementing the ALI method 
in the transfer of the lines which are the most 
difficult to converge. 



Table 1 

Measured Equivalent Widths (reflected spectra) 



Ion 


hu(cV) 


EW''(oV) 


EW''(oV) 


Fc XXVI Lya 


6957.7 


28.06 


17.15 


Fe XXV 


6667.7 


237.7 


145.6 


Fe K (y ^ XIV) 


6400.0 


6.362 


180.6 


S XVI 


3342.8 


4.594 


6.056 


S XVI Lya 


2611.3 


42.90 


47.16 


S XV 


2450.3 


32.28 


42.51 


Si XIV Lya 


1999.8 


27.18 


30.67 


Si XII 


1865.0 


4.892 


7.637 


Fc XXIV 


1673.2 


27.44 


21.00 


Fc XXIV 


1495.6 


16.01 


10.57 


Mg XI 


1343.3 


6.034 


11.19 


Fc XXII-XXIII 


1127.1 


3.698 


6.986 


Fc XXIV 


1110.0 


35.03 


26.44 


No X Lya 


1020.3 


13.46 


16.30 


Fc XVII 


821.10 


4.957 


4.025 


Fc VII 


729.33 


1.951 


2.978 


O VIII Lya 


653.35 


26.37 


32.75 


Fc XXIV 


563.57 


2.466 


2.001 


N VII Lya 


500.14 


3.895 


5.429 


C VI Lya 


367.36 


5.263 


9.009 


Mg X 


194.61 


1.079 


0.468 


O VIII 


163.33 


0.744 


1.360 


Fc XVIII 


127.82 


0.285 


1.578 


O VIII 


120.98 


1.577 


2.706 


Fc XIX 


118.08 


0.711 


2.327 


Fc XXII 


107.81 


0.878 


3.376 



constant density case 



^ constant pressure case. 
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